Multipair contributions to the spin response of nuclear matter 



E. Olsson 1 ' 2 ' 3 and C. J. Pethick 2 ' 3 
1 Department of Astronomy and Space Physics, Uppsala University, Box 515, SE-75120 Uppsala, Sweden 
2 NORDITA, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 
3 Institute for Nuclear Theory, University of Washington, Box 351550, Seattle WA 98195-1550 

We analyse the effect of non-central forces on the magnetic susceptibility of degenerate Fermi 
systems. These include the presence of contributions from transitions to states containing more 
than one quasiparticle-quasihole pair, which cannot be calculated within the framework of Landau 
Fermi-liquid theory, and renormalization of the quasiparticle magnetic moment, as well as explicit 
non-central contributions to the quasiparticle interaction. Consequently, the relationship between 
the Landau parameters and the magnetic susceptibility for Fermi systems with non-central forces 
is considerably more complicated than for systems with central forces. We use sum-rule arguments 
to place a lower bound on the contribution to the static susceptibility coming from transitions to 
multipair states. 
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I. INTRODUCTION 

The spin and spin-isospin responses of nuclear matter 
with various ratios of protons to neutrons are important 
ingredients in many calculations of importance for as- 
trophysical applications. Among these we may mention 
neutral-current processes such as scattering of neutrinos, 
and emission of neutrino-antineutrino pairs, and charged 
current ones such as the emission of neutrinos and an- 
tineutrinos by variants of the Urea process. For a general 
review of neutrino processes in dense matter, we refer to 
Rcf. §. 

In calculating rates of neutrino emission, absorption 
and scattering processes in dense matter, one of the quan- 
tities of interest is the effective interaction between nu- 
cleons, because this is one of the quantities which de- 
termines how matrix elements of the weak current are 
modified by the medium (see, e.g., Rcfs. One ap- 

proach to determining these is to calculate effective in- 
teractions directly as, for example, in Ref. Q. Another 
approach, which has been very successful for liquid 3 He 
is to make use of results for the long-wavelength ther- 
modynamic properties of the system. To be specific, in a 
system with central interactions, some information about 
the effective interaction may be obtained rather directly 
from the spin susceptibility. For the purpose of illustra- 
tion, let us consider such a system with a single species, 
having spin 1/2. The static, long-wavelength spin sus- 
ceptibility is given by 

* = itF' (i) 

1 + Gro 

where [Iq is the magnetic moment of a particle in free 
space, Go is the isotropic part of the Landau quasipar- 
ticle interaction in the spin-triplet channel, and N(0) = 
m*pF/TT 2 h 3 is the density of states per unit volume at 
the Fermi surface, pf being the Fermi momentum and 
m* the effective mass of a quasiparticle. Thus, knowing 



the magnetic susceptibility and the quasiparticle effective 
mass, one can deduce the Landau parameter Gq. 

Landau's theory of normal Fermi liquids provides a 
framework for describing the long-wavelength response 
of the system ||. In the usual account of the theory, 
a key role is played by conservation laws. For example, 
for interactions which conserve the number of particles, 
addition of a single quasiparticle to the system involves 
the addition of a single particle. The total number of 
particles in the system is therefore given by the total 
number of quasiparticles. For interactions which conserve 
the total spin, the total spin may likewise be calculated 
in terms of the quasiparticle distribution, since the spin 
carried by a quasiparticle is the same as that carried by 
a particle. 

The consequences of conservation laws for matrix el- 
ements of the density operator have been described by 
Pines and Nozieres ||, and Leggett has exploited con- 
servation laws to provide a general description of the 
long-wavelength response in terms of contributions due 
to excitation of single quasiparticle-quasihole pairs, and 
those due to multipair excitations [JtJ . A derivation of the 
results from Landau theory is given in Ref. 

The assumption that interparticle interactions con- 
serve the total spin is an excellent approximation 
for atomic interactions, since the spin-orbit and mag- 
netic dipole-dipole interactions between electronic spins, 
which do not conserve total spin, are generally small. The 
situation is quite different for nuclear forces, because of 
the importance of the tensor interaction due to exchange 
of pions and, to a lesser extent, of rho mesons, and the 
significant spin-orbit contributions to inter-nucleon in- 
teractions. For recent analyses of the nuclcon-nuclcon 
interaction, see Refs. [||. 

The purpose of this paper is to investigate the conse- 
quences of the lack of conservation of total spin for the 
spin response of nuclear matter at long wavelengths. One 
of them is that the spin carried by a quasiparticle is dif- 
ferent from that carried by a bare nucleon. Consequently 
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the coupling of the quasiparticle spin to an external field 
is altered. Such effects are well known for nuclear mag- 
netic moments and for matrix elements of the axial cur- 
rent, and they are reviewed in, e. g. , Ref. jlfl]. They are 
also discussed in the framework of Landau Fermi-liquid 
theory in Ref. @. A related effect is that the quasipar- 
ticle magnetic moment is a tensor, not a scalar. A second 
well-known effect is that the interaction between quasi- 
particles has non-central contributions, and the most fa- 
miliar of these is that due to the tensor force. This has 
been considered in Ref. (l^] and more recent works. A 
third consequence of non-central forces is that at long 
wavelengths there are contributions to the response from 
states with more than one quasiparticle-quasihole pair. 
This latter effect appears not to have been widely ap- 
preciated in work on bulk nuclear matter. In this paper 
we examine the general structure of the spin-density re- 
sponse function and deduce the qualitative behavior of 
matrix elements at long wavelengths. We also indicate 
how the standard Landau-theory expression (|l|) must be 
modified when non-central forces are present. Finally, 
we derive an upper bound on the magnitude of the mul- 
tipair contribution to the long-wavelength susceptibility 
from sum-rule arguments. 

This paper is organized as follows. In Section |l| we give 
a brief introduction to response functions and describe 
quantities relevant for our later considerations. Section 
III gives a discussion of the long-wavelength behavior of 
matrix elements, and how this depends on whether or 
not the corresponding operator obeys a local conservation 
law. The consequences of non-central forces for Landau's 
theory of a normal Fermi liquid is examined in Section 
fv|. In Section [v|, we derive bounds on the contribution 
from multipair states to the static, long- wavelength spin 
susceptibility and use microscopic calculations to esti- 
mate how larg e these are. Concluding remarks are given 
in Section [v|. 



II. RESPONSE FUNCTIONS AND SUM RULES 

If a system initially in its ground state is subjected to 
a perturbation 



H' = OlC a e 



Hermitean conjugate, 



(2) 



the Fourier transform of the linear response of the expec- 
tation value of the operator O q at frequency uj is given 

by © 

< O q > w =x(<z,w)C q) (3) 

where the linear response function is defined by 



X 



(ff.«) = £l°£&; 



(4) 



Here j denotes an excited state and the ground state, 
and tUjo = Ej — Eq, where Ej is the energy of the excited 



state, and Eq that of the ground state. The wave vector 
of the perturbation is denoted by q/7i. From Eq. (f|) it 
follows that the static response function is given by 



uj j0 



(5) 



Another important quantity is the static structure factor 

s ^ = lJ2\°l\%- (6) 



Here N is the total number of particles. A third quantity 
we shall find useful, because it can be calculated directly 
from knowledge of the ground-state wave function, is the 
frequency-weighted sum 



W(q) 



U)j 
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|0q[ff,0 q ]|O), (7) 



where H is the Hamiltonian in the absence of the per- 
turbation given in Eq.(||). As we shall demonstrate later, 
the static structure factor and frequency- weighted sum at 
long-wavelengths are useful diagnostics for contributions 
from multipair states. 



III. CONSERVATION LAWS AND 
NON-CENTRAL FORCES 

We now investigate the effect of non-central compo- 
nents of the nuclear force on the linear response of the 
system. One of the reasons for the spectacular success 
of Landau Fermi- liquid theory in providing a framework 
for the quantitative description of long-wavelength low- 
frequency properties of liquid 3 He at low temperatures is 
that the physical observables of greatest interest exper- 
imentally correspond to quantities such as the particle 
density, the spin density, and the particle current den- 
sity, which satisfy local conservation laws. To make this 
point explicit, let us imagine that the operator O(r) sat- 
isfies the conservation law 



dO(r) 
dt 



V-j(r) =0, 



(8) 



where j(r) is the associated current density. On taking 
the matrix element of the Hermitian conjugate of this 
equation between an excited state j and the ground state 
and Fourier transforming in space, one finds 



(9) 



This shows that the matrix element of the operator sat- 
isfies the condition 



q ■ (jq)jO 

LU j0 



(10) 
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Thus for a state for which uijo tends to a nonzero value 
as q tends to zero, the matrix element of the operator 
tends to zero in this limit, provided only that the matrix 
element of the current remains finite. 

The expression ( |To| ) provides the basis for an investiga- 
tion of contributions to physical quantities coming from 
single-pair states and multipair ones. As an example, let 
us consider the case of density fluctuations, that is O q 
is equal to the Fourier transform of the particle density 
operator p q . For a single-pair state, with a quasihole in 
the state with momentum p + q and a quasiparticle of 
momentum p, the matrix element of the current operator 
at long wavelengths is proportional to p/m, where m is 
the mass of the particle. Since the excitation energy of a 
single pair is uj j0 ~ e P + q - e p = v p .q, where v p = V p e p , 
the matrix element of the density operator tends to a 
constant. 

Let us now examine the contributions of single-pair 
states to the static density-density response function, 
the static structure factor, and the frequency-weighted 
sum. In a degenerate Fermi system at zero temperature, 
the number of states with a single quasiparticle-quasihole 
pair having total momentum q is of order nq/pp per unit 
volume for g < pp. The typical energy of a single pair 
excitation is of order VFq, where vf is the Fermi velocity, 
and consequently the contribution to the static structure 
factor varies as q, that to the static response function as 
q°, and the contribution to the frequency- weighted sum 
as q 2 . These behaviors are the same as for a free Fermi 
gas. Similar results apply for the corresponding response 
function for the spin density, irrespective of the nature 
of the interaction. We shall not discuss possible collec- 
tive modes here, but the matrix elements of the density 
operator to such states, which have energies sq, where s 
is the velocity of the collective mode, may be shown to 
vary as q 1 ! 2 . Thus the contributions to x(<?7 0), S(q) and 
W(q) have the same dependence on q as those from the 
single-pair states. 

Now let us turn to the multipair contributions. As 
Pines and Nozieres demonstrated p^j , from Eq. ( ^p| ) it 
follows that the matrix element of the density operator 
vanishes at least as rapidly as q for q — ► 0. If the system is 
translationally invariant, the ground state of the system 
is an eigenstate of the total current, and from this it fol- 
lows that matrix elements of the current between states 
which are orthogonal in the limit q — > must vanish in 
that limit. If (j q )jo can be expanded in a Taylor series, 
this implies that (pq)jo varies as q 2 for small q. The 
density of multipair states does not depend crucially on 
q for small q, and therefore it follows that the contribu- 
tions of multipair excitations to x(g, 0), S(q) and W(q) 
vary as q A for translationally invariant systems, and as 
q 2 for systems without translational invariance. 

Now let us turn to the spin response. If the interaction 
conserves total spin, the contributions to the moments of 
the spin-spin response function have the same behavior 
as for the density. The multipair contributions vary as q 2 , 
since the spin current is not generally conserved, except 



for very special forms of the interaction. For nucleons, 
the interaction has non-central components and therefore 
spin is not conserved. Consequently, matrix elements 
of the spin density operator do not vanish in the limit 
q — » 0, and one arrives at the striking conclusion that the 
contributions to x(9j0), S(q) and W{q) are all nonzero 
for q — > 0. In the following section we shall examine 
the physical content of this result in terms of Landau's 
theory of normal Fermi liquids. 

IV. LANDAU THEORY 

In this section we discuss the long-wavelength low- 
frequency response of the system within the framework of 
Landau Fermi-liquid theory. The basic ideas have been 
presented in the paper by Leggett 0, and we shall ex- 
tend his treatment to the case of tensor forces. While 
Leggett 's work was couched in terms of field-theoretical 
methods, we shall follow more closely the phenomenolog- 
ical approach of Landau. 

Haensel and Dabrowski have shown that for a magnet- 
ically polarized Fermi liquid, the tensor force causes the 
Fermi surface to become aspherical 0] . In this paper we 
shall confine our attention to small polarizations, so the 
distortions of the Fermi surface may be treated as small. 

In a magnetically polarized Fermi system the quasi- 
particle distribution function is a matrix (n p )a/3 in spin 
space. The quasiparticle energy, which is likewise a ma- 
trix in spin space, is defined as the functional derivative 
of the total energy with respect to the quasiparticle dis- 
tribution: 

8E[n p ] = ^(e p ) a /3(Sn p )p a . (11) 

a/3 

The quasiparticle interaction is a matrix in two pairs 
of spin indices, and it is given as the second functional 
derivative of the energy with respect to the distribution 
function 

S 2 E[n p ] = g (/p,P')a/3 : a'/3'(<5n. p W(<H>')/3'a'- (12) 

The magnetic moment of a quasiparticle is related to 
the dependence of the quasiparticle energy on the mag- 
netic field 7i. It is convenient to write the quasiparticle 
energy in the form 

3 

(e p )a/3 =e P ^a/3+X) e P <7 ^' ( 13 ) 

»=1 

where the a 1 are the Pauli matrices. When non-central 
forces are present, application of a magnetic field will 
tend to produce a polarization of the quasiparticle distri- 
bution for a particular momentum in a direction different 
from that of the applied field. Thus the magnetic moment 
of a quasiparticle is a tensor, and it is given by 
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(14) 



If only central forces act, and one neglects excitation of 
degrees of freedom other than the nucleons, the total spin 
is conserved, and the magnetic moment is proportional 
to the component of the spin in the direction i. It then 
follows from conservation of the component of the total 
spin in the direction of the magnetic field, that the mag- 
netic moment of a quasiparticle in a single-component 
Fermi system is equal to the moment of a bare particle. 

The form of the magnetic moment tensor when non- 
central forces are present may be written down imme- 
diately from the requirement that it be invariant under 
simultaneous rotations of the directions of the magnetic 
field and of the momentum of the quasiparticle. This 
dictates that the expression for the magnetic moment is 



fiij(p) = My + 



2 



PiPj 

P 



(15) 



Here p., which is not generally equal to the bare moment 
/Xo, is the magnetic moment averaged over directions of 
the quasiparticle momentum and p,j> is a coefficient char- 
acterizing the strength of the off-diagonal part of the 
magnetic moment tensor. For the specific case of a one- 
pion-exchange potential this form was derived long ago 
by Miyazawa |u| . The general form of the quasiparticle 
interaction may be determined by symmetry considera- 
tions. For purely central interactions, and for an unpo- 
larized system the interaction must be invariant under 
rotation of the spins of the two quasiparticles, and there- 
fore the interaction must be of the standard exchange 
form 



f per,p' < 



/p 



ffpp' 



(16) 



where we have used a compact matrix notation as usual. 
With the matrix indices included, this equation is 

(/p CT ,p'<T')a/3,a'/3' = /pp'<W<W' +5pp'(°')a/3 " (cr') a 'p>. 

(17) 

For calculating low-temperature properties, the quasipar- 
ticle interaction is needed only for p and p' equal to the 
Fermi momentum p-p , and in this case / and g are func- 
tions only of the angle between p and p', which we denote 
by 9. 

When non-central forces are present, the interaction is 
not invariant under rotations of the spins and the mo- 
menta separately, but only under simultaneous rotations 
of spins and momenta. Then there can be additional 
terms of the form erp er'-p, and other terms with one or 
both of the momenta p in this expression replaced by p'. 
For nuclear matter, the most important contribution to 
the non-central part of the interaction comes from pion 
exchange, so it is generally assumed that the interaction 
may be taken to be of the form suggested by the result 
one finds for one-pion exchange 



/p<T, P ' CT ' = / P p' + g pp >crcr' + fcp p /(3cr u er'-u - er-er'), 

(18) 

where u = p — p'. As for / pp ' and g pp i the function k pp i 
is also a function only of 9 at low temperatures. In earlier 
work, e.g. Refs. |l2||lE| ], the non-central contribution to 
the quasiparticle interaction has usually been expressed 
in terms of a function h which is related to k by the 
equation k = hu 2 jp\. As we shall explain elsewhere, the 
advantage of the parametrization we use is that the ex- 
pansion of k in terms of Legendre polynomials converges 
much more rapidly than that for h |l(|. Therefore, all 
of these functions can be expanded in terms of Legendre 
polynomials of cos 9. For example, for k we write 



k pp' = ^kPiicost 



(19) 



1=0 



Now we quote results for the static magnetic suscep- 
tibility. Previously, Haensel and Jerzak jl7| performed 
a similar calculation including the tensor contribution 
to the quasiparticle interaction but neglecting multipair 
contributions and the renormalization of the magnetic 
moment. The magnetization of the matter is 



Mi 



p 



Tr p,i j a jSn p + AI\ 



(20) 



where Mm is the contribution to the magnetization com- 
ing from multipair excitations. And the spin susceptibil- 
ity is defined as 



X 



(21) 



H=0 



For simplicity, we here give the result obtained if one in- 
cludes only Landau parameters with I < 2 for the central 
part of the interaction, and takes only the I = term 
in the tensor part of the interaction k. Using Eq. ( p^| ) 
for the magnetic moment and including the tensor in- 
teraction up to second order, we arrive at the following 
expression for the static spin susceptibility: 



p 2 N(0) 



1 + Go - K 2 



HtirN(p)K 
1 + Go 



- N(0)» T 2 



(22) 



The corresponding result if renormalization of the mag- 
netic moment and multipair contributions are neglected 
is x = ^lN(0)/(l + G - K 2 /8). More details about 
the above calculation and the choice of parametrization 
will be given in a later paper, fig] . Next we shall obtain 
a lower bound on the magnitude of the last term in Eq. 
(^2|), the contribution from multipair excitations. 
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V. A BOUND ON MULTIPAIR CONTRIBUTIONS 
TO THE SPIN SUSCEPTIBILITY 

In this section we shall consider the part of the static 
susceptibility that cannot be calculated within Landau 
theory, the contributions from multipair excitations. Our 
approach will be to employ sum-rule arguments to place 
a lower bound on the quantity. Since the excitation en- 
ergies of the system are positive, it follows for any choice 
of the energy fl that 



El <^ql°> 



2 (0 - U)j ) 

UJjQ 



> 0. 



(23) 



long-range correlations. The difficulty may be illustrated 
by considering experience with calculations of the static 
structure factor for density response for liquid 4 He. On 
the basis of general arguments of the sort given above, the 
static structure factor should vanish in the limit q — > 0. 
However, early calculations gave structure factors which 
were nonzero in this limit. Our conclusion is that sum 
rule arguments of the sort given above can in principle 
provide valuable information about multipair states, but 
that more detailed calculations of structure factors are 
needed before they can provide useful quantitative esti- 
mates. 



If one takes f2 to be the mean excitation energy uj 
W(q)/S(q), this yields the inequality 



2nS(q) 



(24) 



which is valid for all q. Let us now appl y th is result in 
the limit q — > 0. As we argued in Sec. Ill, single-pair 
excitations and collective modes do not contribute to x 
in this limit, and therefore it follows that 



Xm(<7,0) > 



2nS(q) 



(9^0), 



(25) 



where xm(<Z, 0) is the multipair part of the static response 
defined in Eq. (||). 

We now apply this result to pure neutrons, using values 
of the static structure factor and the mean excitation 
energy in the limit q — * taken from Ref. jl8). In that 
reference the authors calculated the spin-spin response 
function, while in this paper we have worked in terms of 
the response of the magnetization density. The response 
functions we need are therefore /ip times those of Ref. 
& 

At the saturation density of nuclear matter, n = 
0.16 fnT 3 , these are Sm(?=0) ~ 0.19/ig and u(q=0) w 
63 MeV. For the total susceptibility we take the value 
X ~ 0.38xf from the recent calculations of Fantoni et 
al. using the auxiliary field diffusion Monte Carlo 
method. Here xf = 3n^Q/2ep is the susceptibility of a 
free Fermi gas. The quantity £f is the Fermi energy of 
the non-interacting gas, and at a density n = 0.16 fm -3 
this is approximately 59 MeV for neutrons. The value 
for the susceptibility obtained in Ref. |l!| is close to the 
results of earlier calculations which included fewer corre- 
lations, (See e.g. Ref. Thus from Eq. @ we find 
Xm(<Z, 0)/x > 0.63. If this result were to be taken at face 
value, it would indicate that multipair states contribute 
more than ~ 60% of the total static response function. 
This would appear to be unrealistically high. The calcu- 
lations of Ref. |l8| were designed primarily to investigate 
the spin response at finite wavelengths, at which phe- 
nomena such as pion condensation would be expected 
to occur. Estimating long-wavelength response is a diffi- 
cult problem because it requires a careful investigation of 



VI. CONCLUSIONS 

One of the main conclusions of this paper is that, when 
non-central forces are present, there are contributions 
to the static magnetic susceptibility at long wavelengths 
that cannot be calculated within the framework of Lan- 
dau Fermi-liquid theory. In addition, non-central forces 
renormalize the magnetic moment of a quasiparticle, and 
introduce additional terms in the expression for the effec- 
tive interaction between quasiparticles. Consequently, it 
is much less straightforward to obtain information about 
Landau parameters from the magnetic susceptibility than 
it is when there are only central forces. This may be seen 
by comparing Eqs. (|l]) and ([22]). 

Among quantities which need to be understood better 
in order to relate quasiparticle interactions to the mag- 
netic susceptibility are the strength of transitions to mul- 
tipair states and the renormalization of the magnetic mo- 
ment of a quasiparticle. In order to place more reliable 
bounds on the strength of transitions to multipair exci- 
tations it would be valuable to have improved estimates 
of the static spin structure factor at long wavelengths. 
The renormalization of the magnetic moment has been 
studied recently by Cowell and Pandharipandc using the 
correlated basis approach developed from the variational 
methods used for quantum liquids pl| . They find that 
the magnitudes of magnetic moments are reduced by ap- 
proximately 10% for nuclear matter with proton fractions 
ranging from 0.2 - 0.5 and for densities between 0.08 fm -3 
and 0.24 fm" 3 . 

In this paper we have focussed mainly on a single- 
component system of fermions with spin- 1/2. Our ar- 
guments may be extended straightforwardly to multi- 
component systems, such as nuclear matter with an ar- 
bitrary ratio of neutrons to protons. In addition, similar 
arguments may be applied for the spin-isospin response. 
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